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Acoustic  Backscattering  From  the  Basin  and 
Margins  of  the  Arctic  Ocean 
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Sound  waves  created  by  high-energy  explosives  were  used  to  measure  teverberation  and  hackscat¬ 
tering  in  the  Arctic  Ocean.  From  an  ice  camp  in  the  Beaufort  Sea.  signals  were  hackscatlcred  from  the 
continental  margin  and  other  major  features  of  the  basin.  An  acoustic  array  was  used  to  analyze  the 


signals  in  azimuth.  Based  on  these  data  we  constructed  charts  of  normalized  backscatter  level,  which 
can  be  compared  with  known  topography  of  the  Arctic  Ocean.  Resolution  of  this  remote  sensing 
technique  for  our  experimental  conditions  (analysis  frequency.  9  Hz;  averaging  time,  20  s;  array  size. 
600  x  600  m)  is  about  8?  in  azimuth,  15  km  in  radial  extent,  and  I  km  in  depth.  We  obtained 
interpretable  signals  out  to  2700  km  For  1000  km  or  less,  the  backscatter  charts  provide  continuous 
coverage  of  the  margins  and  major  submerged  features.  Strong  returns  from  about  73.2rN.  179. (TW 
indicate  one  or  more  prominent  features  that  many  contemporary  topographic  charts  do  not  show .  but 
which  we  believe  to  be  real.  Other  strong  returns  are  evident  from  the  Northwind  Escarpment  and 
from  the  continental  slopes  of  Alaska  and  the  Canadian  Archipelago.  Via  a  backscatter  model,  wc 
estimate  the  roughness  pioduct  (i  nis  height  times  correlation  radius)  of  these  features  to  be  about  4500 
m2  on  average,  f  |  . 


1.  Introduction 

In  March  and  April  1978,  an  ice  camp  was  established  and 
occupied  in  the  Arctic  Ocean  to  remotely  measure  several 
topographic  and  geophysical  properties  ot  the  ocean  bound¬ 
aries.  The  camp  was  located  approximately  at  73.5°N  and 
150. 5°W  in  the  Beaufort  Sea  over  the  Canada  Abyssal  Plain, 
where  the  water  depth  is  about  3800  m.  We  report  here  on 
acoustic  backscattering  from  the  entire  ocean  basin  and  its 
margins,  the  major  features  of  which  are  shown  in  Figure  I. 

The  experiment,  known  as  CANBARX  (Canada  Basin 
Acoustic  Reverberation  Experiment),  used  a  two-dimen¬ 
sional  horizontal  array  to  receive  the  backscattered  acoustic 
signals.  Sound  was  created  at  the  camp's  location  with  high- 
energy  explosives  and  took  as  long  as  3600  s  to  reach  and 
return  from  physiographic  features  2700  km  away.  Interest 
in  such  backscattered  signals  is  two-fold.  First,  the  ocean 
can  be  used  as  a  long-range  acoustic  communication  chan¬ 
nel,  and  the  effect  of  features  such  as  shown  in  Figure  I 
needs  to  be  understood.  Second,  the  backscattered  signals 
carry  information  on  the  nature  and  location  of  these  fea¬ 
tures,  thus  affording  wide-area  knowledge  from  a  single 
measurement  location,  including  the  possibility  of  identify¬ 
ing  unknown  features  and/or  features  of  unusual  relief.  As 
will  be  seen,  the  CANBARX  array  (and  our  processing 
method)  averages  over  or  ignores  certain  features,  but  in 
demonstrating  the  attributes  of  the  experimental  method  we 
hope  to  point  the  way  to  future  higher-resolution  experi¬ 
ments,  as  well  as  to  present  some  new  physical  results. 

Acoustic  backscattering  (reverberation)  on  a  basin-wide 
scale  has  been  previously  measured.  Mcllen  and  Marsh 
f  1 963]  noted  that  ‘reverberation  in  the  Arctic  is  observable 
for  no  less  than  30  minutes  after  an  explosive  detonation  of 


1  Also  at  Department  of  Electrical  Engineering  and  Computer 
Science,  Massachusetts  Institute  of  Technology,  Cambridge,  Mas¬ 
sachusetts  02139. 

1  Also  at  Department  of  Ocean  Engineering,  Woods  Hole  Ocean¬ 
ographic  Institution.  Woods  Hole,  Massachusetts  02543. 

1  Now  at  Applied  Physics  Laboratory,  1  he  Johns  ilopkms  Uni¬ 
versity,  Laurel.  Maryland  20810. 

Copyright  1982  by  the  American  Geophysical  Union. 


moderate  size."  Kutschale  [1969]  observed  the  effects  of 
bottom  topography,  also  in  the  Arctic,  on  both  propagation 
and  reverberation.  Reverberation  in  the  North  Atlantic 
Basin  was  measured  and  modeled  by  Gocnner  and  BUitstein 
[1975]  and  Blalstein  [1978],  again  with  the  use  of  explosive 
sources.  The  essential  distinction  between  the  previous  and 
present  work  is  our  use  of  an  array  to  analyze  signals 
received  as  a  function  of  horizontal  angle.  Such  spatial 
analyses  support  identification  of  specific  features  in  both 
azimuth  and  range  and  provide  detail  heretofore  unavailable. 

In  the  course  of  CANBARX,  other  important  Arctic 
properties  were  measured,  including  ambient  noise  and,  via 
refraction  and  reflection  experiments,  bottom  phase  speeds 
and  layer  thicknesses.  Papers  are  in  preparation  and  in  press 
on  these  results  [Baggeroer  and  Falconer,  1982]. 

2.  Data  Acquisition  and  Processing 

The  camp  was  established  on  new  ice  approximately  2  m 
thick,  with  pressure  ridges  and  old  ice,  and  hence  much 
greater  thicknesses,  about  1-3  km  away,  depending  upon 
direction.  Holes  were  drilled  in  the  2  m  ice  for  the  array, 
which  was  arranged  as  shown  in  Figure  2.  All  1 1  hydro¬ 
phones  of  the  array  were  suspended  to  a  depth  of  61  m  and 
were  spaced  logarithmically  in  the  horizontal  along  two 
orthogonal  legs.  The  array’s  aperture  is  somewhat  over  five 
wavelengths  for  data  at  9  Hz  to  be  presented  here. 

Recordings  were  made  for  up  to  4800  s  after  the  shot 
instant  on  a  12  channel  digital  acquisition  system;  a  typical 
channel  is  schematically  illustrated  in  Figure  3.  An  important 
feature  is  the  gain  ranging  amplifier  which  provided  a 
dynamic  range  of  120  dB,  enabling  acquisition  of  very  low 
level  signals  backscattered  from  remote  areas,  as  well  as 
higher  levels  from  closer  regions.  Through  its  rapid  re¬ 
sponse,  it  also  enabled  acquisition  of  lower  levels  in  the 
vicinity  of  a  strong  scatterer,  along  w  ith  the  high  level  from 
the  scatterer. 

Some  processing  was  done  at  the  camp  to  assure  adequate 
data  quality  and  to  help  guide  the  experimental  program.  The 
bulk  of  the  data,  however,  was  processed  thereafter.  Proc¬ 
essing  steps  are  shown  schematically  in  Figure  4.  For  results 
presented  here,  the  data  were  analyzed  in  3  Hz  bands 
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Fig.  I.  Major  features  of  the  Arctic  Ocean  basin.  Some  of  them,  with  approximate  coordinates,  are  Barrow 
(Alaskan)  coast  (75°N.  I40°-I60'W).  Canadian  archipelago  (78°-83°N.  60°-l25'W|.  Greenland  coast  (83°N,  )0°-60°W), 
European  and  Asian  margin  (80°N,  10°  I60°F.),  Chuckchi  Cap  (78°N.  I600~I80°W),  Northwind  Escarpment  (75°-79°N. 
I55°W),  Alpha  ridge  (85°N.  90°-l80°W).  Lomonosov  ridge  (along  about  80°N.  I40°E  through  the  North  Pole  to  about 
85°N,  60°W).  Redrawn  from  Canadian  Hydrographic  Service  [  1979).  The  encircled  cross  is  the  source/array  location  of 
CANBARX. 


centered  on  9  Hz,  and  beamforming  was  accomplished  via 
the  maximum  likelihood  method  [Baggeroer,  1978;  Capon. 
1979;  McDonough,  1979]  to  yield  an  azimuthal  resolution  of 
about  8°  (0.14  rad)  for  the  signal/noise  ratios  typically 
encountered. 

The  beamformed  sound  signals  are  highly  time  dependent, 
due  not  only  to  the  impulsive  sound  source  but  also  to  highly 
variable  hackscattering  in  time  (range).  With  a  signal  sam¬ 
pling  rate  of  ?<'0  Hz,  we  formed  estimates  of  the  time- 
dependent  pressure  spectral  density  by  averaging  the  signal 
over  both  the  3  Hz  bandwidth  and  a  20  s  time  window.  The 


time  window  was  moved  forward  in  4  s  steps  and  gave  a  set 
of  data  points  versus  time  for  each  bearing.  Finally,  the 
bearing  was  swept  through  360°  in  6°  steps,  to  yield  an  array 
of  data  points  pertinent  to  basin  backscattering.  Reasons  for 
and  implications  of  the  choice  of  the  20  s  time  window  and 
the  9  Hz  center  frequency  will  be  explicitly  discussed  in  this 
paper. 

To  genciaic  impulsive  sources  of  sound,  we  exploded 
TNT  at  a  pressure  depth  of  244  m,  w  ith  a  range  in  mass  from 
50  to  400  kg,  mostly  at  200  kg.  A  source  level  monitors 
occupied  the  12th  channel  of  the  data  acquisition  system. 
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using  a  hydrophone  of  much  lower  sensitivity  than  the  other 
II  in  the  array.  But  the  monitor  hydrophone  was  deployed 
improperly  just  below  the  ice  making  interpretation  next  to 
impossible.  It  was  too  close  to  a  pressure  release  surface  and 
a  region  of  layered  cavitation  induced  by  the  intense  reflect¬ 
ed  signal  [Malme  cl  al.,  1964;  Gaspin  a  uL,  1979] .  Instead  of 
measured  source  data,  we  rely  therefore  upon  standard 
values  for  underwater  TNT  explosions,  the  most  compre¬ 
hensive  results  coming  from  the  formulas  and  charts  of 
Weston  [I960]  and  Urick  [1975],  Values  taken  from  these 
publications,  and  detailed  intercomparisons  of  CANBARX 
backscatter  data,  indicate  that  the  maximum  uncertainty  in 
source  level,  with  the  use  of  this  procedure,  is  ±3  dB, 
including  all  other  uncertainties  in  signal  level  (due  to 
calibration,  voltage  standards,  etc.).  The  relative  uncertainty 
in  our  data,  however,  is  less  than  I  dB. 

In  many  ways  the  Beaufort  Sea  ice  cover  in  March  and 
April  is  an  excellent  platform  for  scientific  work.  It,  howev¬ 
er.  is  known  to  drift  under  action  of  wind  and/or  current 
stress.  During  our  experiment  the  camp  remained  within  an 
ellipse  of  about  4  km  N-S  semi-axis,  and  7  km  E-W  semi¬ 
axis.  It  also  rotated  with  a  standard  deviation  of  about  2°. 
Both  motions  appear  random.  We  have  not  corrected  our 
backscatter  data  for  these  motions,  since  they  fail  well 
within  other  uncertainties  in  or  averages  over  the  data,  and 
in  any  case  have  small  effect  in  relation  to  the  wide-area 
perspective  of  our  experiment. 

3.  Reverberation  Data 

We  form  the  reverberation  level 


RLt  =  10  log 


T„(S(f0,  1)) 
A(/0) 


dB  re  1  m 


(I) 


where  47rA(/0)/pc  is  equal  to  source  energy  (p  water  density 
and  c  sound  speed);  5(/0,  t)  is  the  spectral  density  of  p2(t), 
the  sound  pressure  squared;  Tp  is  the  processor  averaging 
time;  and  where  {  ■  •  •  )  indicates  the  average  over  both  Tp 
and  bandwidth  Wp. 

In  our  study  we  chose  values  of  f0  =  9  Hz,  Wp  =  3  Hz, 
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Fig.  3.  Recording  system  for  each  channel 

;  20  s.  To  the  exte 
distributed  over  at  least  20  s. 


and  Tp  =  20  s.  To  the  extent  that  S(/0.  /)  is  uniformly 


RL,  .  2o  =  Rl- :o  +  10  log 


20 


(2) 


and  subsequently  we  shall  argue  that  this  is  so,  at  least 
approximately,  for  the  signals  of  interest. 

The  integrated  source  spectral  density  at  the  frequency  /0 


A(/„) 


.9,(/0.  t)  di 


(3) 


where  Ts  is  the  duration  of  the  explosion  and  5,  is  the 
spectral  density  evaluated  at  the  refeience  distance  of  1  m 
from  the  source.  The  explosion  duration  7j  is  about  0.3  s.  but 
as  indicated  earlier  we  use  standard  values  of  A(/u)  rather 
than  perform  the  integration  in  (3)  with  measured  values  of 

S5(fo,  »■ 

A  three-dimensional  view  of  RL2 o  versus  bearing  and 
time.  Figure  5,  shows  significant  structure  in  the  backscat- 
tered  sound  field  out  to  the  maximum  displayed,  3600  s  (3431 
s  after  shot  instant).  The  level  at  the  shot  time  (169  s)  is  not 
physical,  representing  saturation  of  the  sensitive  hydro¬ 
phone-preamplifier  sensors  of  the  array.  Thereafter,  the 
sensors  recover  and  display  RL 20,  initially  with  significant 
structure  through  360°.  and  then  increasingly  with  confine¬ 
ment  toward  the  northerly  sector  as  the  signals  return  from 
more  northerly  regions  of  the  basin.  (Curiously,  Figure  5 
shows  noise,  at  say  180°  true  and  3600  s,  to  be  less  than  that 
before  the  shot,  an  effect  often  observed  and  to  be  described 
more  fully  in  another  paper  on  ambient  noise.  In  brief,  for 
low  ice  stress  and  therefore  low  noise  conditions,  an  explo¬ 
sion  relieves  stress  and  reduces  the  noise  nearby  for  some 
hours  thereafter.) 

Contour  diagrams  of  RL 20  are  more  suitable  for  quantita¬ 
tive  interpretation.  An  example  is  shown  in  Figure  6.  With 
the  aid  of  the  overlay  of  the  basin's  major  features,  one  can 
discern  significant  returns  from  the  Barrow  Coast  to  the 
south,  from  the  Canadian  Archipelago  to  the  southeast  and 
east,  and  from  the  Northwind  Escarpment  to  the  west.  An 
extensive  set  of  reverberation  diagrams,  from  which  the  one 
shown  here  is  taken,  is  available  [ Zittel ,  1979] 

With  appeal  to  linearity,  a  reverberation  diagram  such  as 
Figure  6  can  be  applied  directly  to  a  wide  class  of  acoustic 
communication  problems  [Urick,  1975],  It  will  become  obvi¬ 
ous  in  the  next  section,  however,  that  such  results  are  in 
general  specific  to  f0.  to  the  sourcc/recciver  depths,  and  to 
their  geographic  location. 


4 


t 

! 


4<i.4i.l'f4'k4«a.l  *1 


\ry.  u  »  w 


Dyer  m  ai  .:  Acoustic  Baikscaihr  From  Ak(.ik  Basin 


CANBARX  Processor 


X.tnAl)  I.lmAO 


Window,  Estimate 

FFJ  - •"  Spec  Jr  oJ  Steering 

and  '  •"  Covononce  ^  Function 


Conversion 
to  PTmer 

a  i_> 

Correc  tion 


p(u/l  *(  d*R_'d  ) 


Fig.  4.  Block  diagram  of  the  array  processor. 
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4.  Conversion  or  Reverberation  to  Backscatti.r 
Data 

Further  processing  of  reverberation  diagrams  is  needed 
when  the  objective  is  better  understanding  of  physiographic 
features.  As  may  be  seen  in  Figure  6,  for  example,  contours 
enclosing  scattercrs  at  a  given  RL^o  niay  extend  over  a 
considerable  range.  Surely,  a  scatterer  at  large  range  is  far 
more  important  than  one  at  small  range  with  the  same  RLX. 
Also,  reverberation  diagrams  give  no  guidance  on  the  depth 
of  important  scatterers.  We  will  show  that  the  CANBARX 
source/receiver  combination  emphasizes  scatterers  in  the 
upper  part  of  the  water  column.  We  will  also  normalize  the 
reverberation  data  to  reduce  range-dependent  effects,  there¬ 
by  facilitating  identification  of  physiographic  features  and 
quantitative  comparison  among  them. 

4.1  One-Way  Propagation  in  an  Ideal  Channel 

The  sound  speed  profile  in  the  Arctic  Ocean  is  approxi¬ 
mately  a  linear  function  of  increasing  depth,  with  some 
deviations  in  the  upper  500  m  or  so  which  we  shall  ignore. 
Thus,  all  acoustic  rays  follow  upward  refracting  paths  that 


are  to  a  good  approximation  circular  arcs.  From  Snell's 
Law,  the  water  depth  Dj  and  source  depth  determine  the 
maximum  ray  angle  0j,  (in  small  angle  approximation): 

,  2 #</>,  -  ;,) 

0i*  „ - —  (4) 

where  g  is  the  sound  speed  gradient,  and  <i  the  sound  speed 
at  the  source  depth.  In  this  model,  all  rays  described  by  |0|  < 
0i  follow  refracted  paths  completely  in  the  water  column, 
reflecting  periodically  from  the  ice  canopy,  while  rays  |0|  > 
ff,  travel  at  least  partially  in  the  bottom.  As  is  customary,  we 
neglect  the  latter  because  bottom  paths  are  known  to  have 
much  higher  attenuations  than  water  paths.  With  g  -  1.85  x 
KT2  s1.  f|  =  1.45  km/s,  and  D\  -  Zi  =  3560  m  as  typical 
values,  we  have  0|  =  0.30,  which  is  small  enough  for  use  of 
small  angle  approximations  whenever  convenient.  Also,  at  9 
Hz  the  array  steered  in  the  horizontal  emphasizes  returns 
from  within  |0|  £  0.4,  further  justifying  the  neglect  of  bottom 
interacting  paths. 

We  neglect  cavitation  which  is  expected  at  the  water-ice 
interface  in  the  vicinity  of  the  source,  and  defer  its  discus- 
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Fig.  5.  Reverberation  level.  RL^U  on  a  relative  basis  versus  bearing  (geographic  true)  and  time.  Charge  weight  2(X) 
kg.  source  depth  244  m,  receiver  depth  61  m.  center  frequency  9  Hz.  A  nominal  sound  speed  of  I  .48  kill's  is  used  to  fix  a 
range  scale  (experiment  I  t). 
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Fig.  6.  Contours  of  RL:n  in  dB  re  I  m.  Charge  weight  50  kg  [A(9  Hz)  =  2  x  10”  (uPa:  m  :  :>  Hz  ')  (experiment  31). 
Noncontour  lines  are  200  m  depths  from  known  topography.  The  North  wind  Escarpment  (at  800  m>  is  also  shown. 


sionforthe  end  of  this  section.  Following  Smith  [1971,  1974] 
we  assume  the  source  to  be  a  delta  function  and  write 


S(f,  t)  = 


4  Mf)eaR 


( ~  (7~i  ~  {R&Kx  Mde 

J8„  ~~  X(6 )  tan  10| 


which  is  a  ray  averaged  expression  for  the  sound  field.  In 
this  equation  A70)  is  the  ray  period  (cycle  distance),  R  is  the 
horizontal  range,  a  is  the  attenuation  coefficient  (based  on 
loss  which  is  assumed  a  linear  function  of  ray  angle),  <f>  is  the 
ray  angle  of  the  6  ray  as  it  travels  through  the  water  column 
and  is  given  by 

<t>2  =  e2  -  2g(z-~-—  (6) 

C| 

where  z  is  the  depth  measured  in  the  positive  sense  from  the 
ocean  top  to  the  bottom  D\.  Since  d>can  not  be  complex,  we 
see  from  (6)  that  the  integral's  lower  limit  is  8m  =  0  for 
observation  above  (z  <  Zi)  and  9m  =  |2#(z  -  Z|)/C|},c  for 
observation  below  (z  >  zi)  the  source. 

We  note  that  (5)  describes  a  pulse  first  arriving  at  a  time 
after  shot  instant  given  by  Tt  -  ROplUy  and  culminating  at 
T |  -  R8„‘l(*'\.  Its  duration  Li  is,  thuefore. 


f T|  is  given  by  (Rlc i)  (I  -  (2gz|/3cl))J.  For  values  of  R  and  z 
of  interest  in  the  CANBARX  data.  Lt  greatly  exceeds  the 
explosion  duration  7j,  thus  justifying  the  delta  function 


assumption.  The  (integrable)  singularities  in  (5)  [Wesion, 
1980]  predict  that  most  of  the  pulse  energy  is  concentrated  at 
its  culmination,  as  experiments  in  the  Arctic  verify  [DiNa- 
poli  el  a!.,  1978].  While  it  might  be  tempting  to  use  actual 
source  shapes  0  in  convolution  with  (5)  to  better 

predict  the  shape  of  S(f,  i)  at  its  peak,  we  need  not  do  so, 
since  ultimately  we  shall  average  S(/,  r)  in  a  way  that 
suppresses  this  detail. 

As  indicated  in  the  foregoing,  cavitation  at  the  top  of  the 
water  column  adversely  affected  our  measurement  of  source 
strength.  The  hydrophone  deployed  for  this  measurement  at 
times  was  in  the  cavity  region,  and  the  pressures  recorded 
were  therefore  unreliable  indicators  of  source  properties.  On 
the  other  hand,  in  (5),  which  describes  the  sound  field,  we 
neglected  cavitation  and  we  now  justify  this.  For  shots  used 
in  our  experiment  (50  to  400  kg  TNT  detonated  at  244  m)  the 
spectral  peaks  arc  in  the  3  Hz  band  centered  at  9  Hz.  These 
peaks  are  dominated  by  the  bubble  energy  and,  unlike  the 
shock  wave  energy,  reflect  from  the  top  of  the  water  column 
as  do  infinitesimal  signals  [Caspin  el  al.,  1979).  Thus  (5)  may 
be  properly  applied  to  our  experiment. 

4.2  Effective  Channel  Depth 

The  time  integral  of  (5)  gives  the  energy  profile  in  z,  a  few 
values  of  which  are  given  in  Table  1.  We  see  a  substantial 
bias  toward  energy  in  the  upper  part  of  the  water  column. 
Also,  on  average,  the  energy  levels  above  500  m  are  about 
equal  to  that  found  at  500  m.  (The  energy  profile  predicted 
by  (5)  is  singular  at  244  m,  the  source  depth,  a  well-known 
failing  of  ray  theory.  Wave  theory  would  give  the  proper 
finite  result  but  it  would  not  sensibly  affect  the  definition  of 
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TABI.E  t.  Acoustic  Energy  Depth  Profile  (Approximate) 


Depth,  m 

Energy  Level, 
dB  re  t.evel 
at  500  m 

Gradient  of 
Energy  Level 

500 

0 

-1  dB/100  m 

800 

-2 

-  1  dB/200  m 

1200 

-4 

-1  dB/300  m 

1900 

-6 

- 1  dB/400  m 

Source  at  244  m. 


effective  channel  depth.)  Thus,  we  can  define  an  effective 
depth  O,  such  that  the  region  0  s  z  ^  Dr  contains  the 
acoustic  energy  important  to  backscattering.  The  half-ener¬ 
gy  depth  is  at  about  1000  m,  and  setting  De  =  1000  m  is  a 
natural,  although  crude,  definition  of  the  effective  channel 
depth.  We  note  that  others  have  made  similar  definitions 
[e.g..  Buck,  1981]. 

A  scattering  feature  below  2000  m.  for  example,  would  be 
expected  to  be  irradiated  by  at  least  6  dB  less  energy  than 
one  at  500  m  or  above,  and  hence  such  deep  features  are  not 
likely  to  be  important  in  backscattering  diagrams.  Thus  we 
should  expect  the  continental  margin,  rather  than  mid-ocean 
features,  to  be  the  principal  locus  of  backscattering. 

We  may  associate  with  De  an  effective  source  ray  6e  such 
that  |0|  £  6e  describes  the  region  of  energy  important  to 
backscattering.  For  De  =  1000  m  we  have  8e  =  0. 14,  and  this 
source  angle  replaces  0,  in  the  integral  of  (5)  for  the  incident 
acoustic  field. 

4.3  Wedge  Model  of  Physiographic  Features 

Signals  described  by  (5)  are  scattered  by  any  physical 
feature  departing  from  the  ideal  constant-gradient  constant- 
depth  channel.  We  consider  only  departures  from  the  latter 
and  restrict  attention  to  features  measurable  in  the  scale  of 
basin  dimensions,  such  as  the  continental  margin. 

We  hypothesize  that  a  major  feature  may  be  modeled  as  a 
wedge  of  angle  p,  with  the  incident  plane  of  sound  at  the 
angle  y  to  the  wedge  (see  Figure  7).  The  ray  angle  <t>  at  the 
wedge  then  has  a  grazing  angle  to  the  wedge  of  <f>  +  t/i,  where 
i/<  is  the  in-plane  wedge  angle  determined  by 

'F  =  tan  <!>  =  tan  p  cos  y  (8) 

The  plane  defined  by  /3  is  in  actuality  the  mean  of  a  rough 
surface,  and  we  concentrate  here  on  suitable  characteriza¬ 
tion  of  this  plane  for  basin  backscattering.  The  angles  p  and 
•yean  be  estimated  at  each  location  from  present  knowledge 
of  the  basin  and  its  features.  But  we  prefer  to  consider  them 
as  stochastic  variables  with  estimated  or  assumed  statistics 
and  density  functions.  Ultimately,  data  from  diversely  locat¬ 
ed  backscatter  experiments  might  be  useful  in  determining  p 
and  y  independently,  but  for  present  purposes  we  adopt  the 
stochastic  view  as  the  more  direct  way  of  gaining  an  appreci¬ 
ation  of  basin  backscattering. 

We  assume  B  =  tan  p  has  a  Rayleigh  probability  density 
function  (pdf),  with  mean  pti  and  variance  <rB2  =  pH2  ((4/7?)  - 
J).  The  assumed  pdf  as  such  is  not  crucial,  but  when 
sampling  over  a  large  number  of  wedges,  we  take  it  as  likely 
that  the  pdf  is  skewed  towards  the  smaller  values  of  B.  We 
also  assume  that  the  pdf  for  y  is  uniform  in  the  interval  ( -  nl 
2,  7r/2)  and  hence  that  G  =  cos  y  has  a  pdf  with  pc,  =  21  tr  and 
oa  =  1/2  -  4/7T2.  Again  the  precise  form  of  the  assumed  pdf 
is  not  crucial,  and  is  applicable  principally  as  an  ensemble 
description  of  the  basin’s  physiographic  features.  With  these 


pdfs  it  readily  can  be  shown  that  the  pdf  for  VF  is  half- 
Gaussian,  with  p^  =  2pn!tr  and  try2  -  p/?2(  I  -  Such 

a  statistical  view  of  Arctic  Ocean  features  is  ofTered  not  as  a 
proven  experimental  fact,  but  rather  as  a  plausible  model 
whose  use  can  help  in  interpreting  backscatter  data. 

In  a  simila.  vein,  the  wedge  angle  <p'  in  the  plane  of  a  w  ave 
front  is  given  by 

W  =  tan  i !)'  =  tan  /3|sin  y|  (9) 

and  has  a  pdf  and  statistics  identical  to  (.'lose  of 'I'. 

Locally,  a  wedge  restricts  the  ray  angle  to 

0  s  <f>  +  i/)£(4t/i  0)|,:  (10) 

The  lower  limit  clearly  restricts  rays  from  incidence  beneath 
the  wedge,  while  the  upper  limit  can  be  shown  to  restrict 
scattering  from  the  wedge  to  first-order  scattering  only. 
These  restrictions,  together  with  (6),  describe  the  variation 
of  <t>  with  depth  ;  on  the  wedge. 

One  immediately  can  see  that  backscattering  from  Arctic 
Ocean  physiographic  features  is  governed  by  quite  small  ray 
angles.  If  we  take  p*  -  3  x  10 "2  ( p *  =  1 .7°,  pp  ~  2.7°)  and 
note  from  (10)  that  the  maximum  values  of  (d>  +  iM:  are 
distributed  as  1 1>  |0|,  then  {<t>  +  </i)rlm  is  less  than  1.3  x  10“ 1 , 
i.e.,  7.6°,  for  9  less  than  9e. 

In  addition  to  p  and  y,  the  transverse  width  of  a  wedge.  /. 
affects  the  backscattering  of  sound.  We  hypothesize  that  this 
is  approximated  by  its  longitudinal  width, 

/  =  cot  p  (II) 

where  z„  is  the  wedge's  height  (see  Figure  7).  Our  view  of  a 
physiographic  feature  is  thus  of  a  string  of  wedges,  each 
wedge  characterized  by  its  slope  tan  p  and  its  orientation  y 
referred  to  the  sound  wave  vector.  The  slope  and  orientation 
have  a  horizontally  isotropic  scale  /.  which  itself  can  be 
considered  a  stochastic  variable  via  p  and  zw.  We  shall  not, 
however,  attempt  to  specify  its  stochastic  variation  explicit¬ 
ly. 

4.4  Overlap  of  the  Sound  Field  With 
Physiographic  Features 

Irradiated  wedges  representing  a  section  of  the  continental 
margin  are  potentially  capable  of  scattering  energy  back 
toward  the  receiving  array.  In  our  processing  system,  back- 
scattering  ultimately  is  observed  through  a  bearing  window 
fp  (0.14  rad)  and  a  time  window  Tp  (20  s).  The  overlap  of  the 
acoustic  area  defined  by  these  windows  with  the  horizontal 
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Fig.  8.  Acoustic  sector  overlaying  the  horizontal  projection  of  a 
physiographic  feature. 

projection  of  a  feature  is  shown  in  Figure  8.  From  this,  we 
can  find  the  ratio  F  of  included  or  irradiated  wedge  area  to 
acoustic  window  area,  and  its  various  values  are  given  in 
Table  2.  The  quantity  A z  is  the  vertical  extent  of  the  feature, 
an  estimate  of  which  follows  from  the  vertical  coverage  of 
the  sound  field  and  from  dimensions  of  the  feature  itself. 
(See  Figure  9.)  In  particular,  the  acoustic  field  has  a  lower 
channel  depth  De  =  1000  m,  so  that  the  continental  slope 
normally  would  not  be  sensed  down  to  its  transition  to  the 
continental  rise.  Also  the  slope  extends  up  to  Ds,  its  transi¬ 
tion  to  the  shelf.  Hence,  A z  -  Df  -  Ds ,  and  with  a  nominal 
value  D ,  =  200  m,  we  estimate  A z  -  800  m. 

In  composing  Table  2  we  have  assumed  the  scale  /  of  (1 1) 
to  be  larger  than  the  scale  of  the  irradiated  wedge,  which 
permits  us  to  think  of  a  single  value  of  'F  or  'F'  for  each 
overlap.  This  assumption  is  not  restrictive,  because  ulti¬ 
mately  we  shall  treat  T  both  approximately  and  as  a  stochas¬ 
tic  variable,  the  latter  implicitly  allowing  ¥  (and  ¥')  to 
change  within  the  scale  of  an  irradiated  wedge.  Incidence 
angles  y0  and  y,  in  Table  2  are  given  by 

2Az 


cTp 

tan  y.  =  - 

2  RSP 

(12) 

TABLE  2. 

Values  of  T,  the  Ratio  of  Irradiated  Wedge  Area  to 
Acoustic  Window  Area 

2Az 

T  < - 

cTp 

,  2 ^ 

T  > - 

cTp 

A  z 

R4T  <  — 

1 

Ty<y0 

2Az 

- -  ■  y  >  ?o 

cTpV 

A  z 

RV  >  — 
ie 

A  z 

R&V’ 

2A  z  „ 

- y  <  Yi 

cTpV 

A z 

- - .  y  >  Yi 

with  typical  values  y0  =  37°  (for  'F  =  3o>  =  6.8  x  10  ')  and 
yi  =  6°  (for  R  =  1000  km). 

We  argue  that  Table  2  is  well  represented  for  CANBARX 
by  the  first  column  only 


1 


r 


A; 

fl'F'  <  — 


Ac 

Kf/F' 


Ac 

RV  >  — 


(13) 


Since  2\zl(iTp)  -  5.1  x  10~2.  values  of  T  applicable  to  the 
first  column  of  Table  2  cover  more  than  2  standard  devi¬ 
ations  (o\j,  =  2.3  x  10“"  for  ny  =  3  x  10~2  and  the  presumed 
half-Gaussian  pdf).  Also  F  drawn  from  the  second  column 
approaches  (13)  anyway.  That  is,  for  /fT  ’  small  (correspond¬ 
ing  to  y  — »  0)  and  for  RV  large  (corresponding  to  y  — ►  90°), 
we  see  that  F  tends  to  be  (13). 

From  (13)  we  can  define  a  crossover  range  R0,  based  on 
the  mean  w 


Ac 

Ro  -  - -  (14) 

Thus  when  R  <  R0  the  acoustic  window  has  high  probability 
of  being  completely  filled  by  the  wedge  string,  while  for  R  > 
R0  it  does  not,  suffering  in  effect  an  acoustic  loss.  We  shall 
use  this  result  in  normalization  of  the  reverberation  data  to 
obtain  interpretable  backscatter  results. 

With  (iy  =  jx  10’2  and  the  foregoing  estimate  of  Ac.  the 
crossover  range  R0  is  about  190  km.  Since  Ac  and  depend 
upon  the  measurement  system,  Rv  is  not  solely  a  basin 
parameter. 

4.5  Two-Way  Propagation  Model 

Equation  (5)  has  a  simple  6  integral  which,  in  accordance 
with  section  4.2,  we  integrate  to  6,  rather  than  6,.  We 
average  this  result  in  c  from  0  to  £>r,  principally  because  we 
do  not  have  a  priori  knowledge  of  depth  of  the  relevant 
scatterers,  with  the  result  that  the  incident  field  can  be 
expressed  as  [Smith.  1971). 


S, 


2 a  ee 
TrRDr 


(15) 


where  r,  is  the  transmission  coefficient  of  (he  incident  field, 
where  AIT,  can  be  interpreted  as  the  effective  rate  of  energy 
propagating  in  the  channel  (aside  from  a  factor  4n/pc).  and 
where  T,  is  the  range  dependent  time  spreading  of  the 
impulse.  The  latter  is  equal  to  about  I  s  per  100  km  for  the 
full  Arctic  acoustic  channel  (A/  from  (7))  but  is  about  I  s  per 
500  km  for  the  more  energetic  components  in  the  depth 
interval  0  <  c  <  De.  Data  supporting  this  estimate  can  be 
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found  in  Kutschale  [1968]  and  Milne  |I967).  Simply,  S,  is  the 
average  spectral  density  in  the  effective  time  interval  7,. 

The  component  of  the  intensity  vector  normal  to  a  scatter- 
er,  modelled  in  macroscale  as  a  wedge,  is  5,  sin  tpUftc).  Thus 
the  power  flowing  normal  to  the  wedge.  W,.  is  S  .S’,  sin  iMpc) 
where  S  is  the  wedge  surface  area  occupied  b>  (he  process¬ 
ing  windows.  The  energy  is  7,W,.  a  certain  fraction  b  of 
which  is  reradiated  as  backscattered  energy  bTrWr  That  is, 
we  take  each  wedge  to  have  roughness  or  facets  which  are 
the  actual  backscatter  mechanisms,  with  j3  defining  the  mean 
scattering  surface.  While  such  backscattering  is  known  to  be 
highly  dependent  upon  angle  between  mean  scattering  sur¬ 
face  and  wave  direction,  we  represent  the  process  simply  by 
the  parameter  b,  with  justification  that  the  wave/wedge 
interaction  naturally  averages  over  the  spread  of  incident 
and  reflection  angles. 

The  rate  of  energy  backscattered  in  the  time  window  Tp  is 


dEh  =  bT,W, 
dt  T„ 


(16) 


where  time  spreading  caused  by  the  wedge's  range-wise 
extent  is  taken  to  be  uniform,  justified  by  the  presumed 
equality  in  incident  energy  with  depth,  and  again  by  the 
inherent  averaging  over  both  incident  and  reflection  angles. 
We  note  that  the  uniform  hypothesis  requires  T,,  <  2Ac/(t  ffO. 
which  holds  in  our  system  for  ¥  <  5  x  10  \  but  even  when 
violated  still  represents  the  time  spreading  with  acceptable 
accuracy.  We  see  that  wedge  time  spreading,  as  measured 
by  Tp.  dominates  channel  spreading  T,  (and  also  source 
duration  7,),  and  justifies  (2).  That  is,  scattering  and  trans¬ 
mission  back  to  the  array  can  be  assumed  with  comfort  to  be 
approximated  in  the  time  domain  as  a  rectangular  function, 
with  smearing  due  to  radiation  and  transmission  occupying 
an  interval  reasonably  small  compared  to  Tp. 

Conversion  of  dEJdt  to  spectral  density  of  the  backscat¬ 
tered  field  Si,  is  straightforward.  Since  Ah  =  pcE^ir  is  the 
equation  connecting  energy  f*  and  strength  A/,  of  a  field 
radiated  over  2 tt  steradians,  we  have  in  analogy  with  (15) 


AbS  sin  ip  2AbS  sin  ip 

VTb  = 


0\ 

2ttTp  '  "  ttTp  \R2D,.D[ 


(17) 


Note  that  5/,  does  not  depend  upon  Tr,  since  the  incident 
spectral  density  and  backscatter  energy  rate  depend  inverse¬ 
ly  and  directly  upon  Tr.  respectively.  It  is  simple  to  show 
that  the  scattering  area  is  given  by 


I  Rjp  cTp 
2  cos  ip  cos  ip' 


(18) 


which  takes  on  various  values  from  P  in  Table  2  or  as 
approximated  for  CANBARX  in  (13). 

With  reference  to  (I)  the  reverberation  signal  r  is  modeled 
as  TpSjA.  with  dimensions  of  inverse-square  length,  and 
from  ( 17)  and  (18)  is 

r  1 


b*  \( RincTpH 


cos  ip' 


A; 


(19) 


where  the  split  values  refer  to  conditions  in  the  same  order 
as  in  (13).  We  see  that  for  small  R  the  reverberation  signal 
varies  as  R  while  for  large  R  it  goes  as  R  Also  at  large 


R  the  signal  does  not  depend  upon  £,,,  the  resolution  being 
A d(Ry').  but  the  resolution  displayed  will  depend  upon  the 
steering  steps  used. 

Pquation  (19)  is  seen  to  justify  (2).  for  it  gives  a  linear 
dependence  of  r  upon  7,,.  Implicitly.  7,,  must  be  no  more 
than  2A;/t(  'l/).  and  we  take  20  s  as  the  largest  practical  value 
for  our  experiments. 

Let  us  write  r  as  the  estimate  of  r  from  measurements  in 
the  time  window  Tp.  Then  the  model  provided  by  ( 19)  can  be 
restated  as 


by 


cos  ip'  R£p  cTr 


( T,  Th ) 


RV 

RoM-'f 


(20) 


in  which  we  have  used  (14)  for  Rt)  and  where  is  interpreted 
as  an  estimate  from  measurements  of  a  backscattering 
quantity.  We  know  all  parameters  on  the  rhs  except  for  ff7'  of 
which  we  do  not  have  a  priori  knowledge.  Thus  we  average 
over  the  angle  parameters  of  the  basin,  with  the  icsult 


Er, 


by 

COS  Ip' 


4  TT 

~R&  TP  (T|  Th) 


1 

RIRo 


R  <  Ro 

R>  Ro 


(21) 


where  p,,  is  the  angle  average  of  the  Ihs  of  (20).  As  a  model  of 
backscattering,  (21)  normalises  the  measurements  to  yield 
physically  comparable  results  throughout  the  basin.  System¬ 
atic  changes  in  p,,  from  some  overall  norm  can  then  be 
ascribed  to  local  ‘abnormalities’  in  b  or  ff7.  although  we  have 
not  as  yet  separated  these  in  our  results. 

4.6  Attenuation  Coefficient 

Parameters  in  (21)  have  been  estimated  in  prior  sections, 
except  for  exp  \  -2aR\  appearing  in  r,  iy,  we  arrive  at  this 
estimate  here.  As  is  customary,  we  write  a  =  4.34a  so  that 
the  loss  coefficient  is  expressed  in  dB  per  unit  distance.  At 
low  frequencies,  loss  mechanism  candidates  are  (I)  boric 
acid/borate  relaxation  [ Fisher  and  Simmons,  1977],  (2)  un¬ 
derice  scattering  [Diachok,  1976],  and  (3)  volume  scattering 
[Mellen  et  al.,  1974], 

At  9  Hz  the  first  is  estimated  to  be  ai  =  6  x  10  6  dB/km, 
with  frequency  dependence  of  Thus  the  first  mechanism 
has  negligible  effect  in  two-way  transmission  out  to  a  range 
of  2500  km.  Uncertainties  associated  with  a\  pertain  to  the 
pH  of  Arctic  waters,  but  are  within  a  factor  of  2  using  other 
oceans  as  a  guide  [ Lovett ,  1980). 

The  second  mechanism  may  well  be  dominant  in  CAN¬ 
BARX,  or  at  least  nonncgligible.  Diachok  modeled  ice 
scattering  as  reflections  from  cylindrical  keels.  With  use  of 
his  model  and  parameter  estimates  (10.5  keels/km.  keel 
depth  4.8  nt,  and  keel  half-width  9.8  m),  we  find  e;  -  3  x 
10  4  dB/km  at  9  Hz.  On  the  other  hand,  Bass  and  Fuks  (1979] 
modeled  loss  in  a  nonrefraeting  parallel-plate  waveguide 
with  Gaussian  roughness  of  rnts  depth  and  correlation 
radius  A,  on  one  surface.  With  use  of  their  formula  adjusted 
to  account  for  refraction,  and  with  estimates  =  2.4  m  and 
A,  =  70  m  ( Wadbams  and  Horne,  1980;  Lawn,'  and  Wad- 
hams,  1979;  Rothrock  and  Thorndike,  1980],  we  get  a;-2  x 
I0~ 5  dB/km.  The  former  estimate  would  have  a  3  dB  effect  at 
2500  km,  the  latter  10  dR,  and  correspondingly  less  at 
smaller  R.  The  frequency  dependence  of  o;  is  /'  (Diachok) 
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Fig.  10.  Baekscattering  levels  for  the  Arctic  basin  (experiment 
14).  Contours:  -52.  -56.  -40  dB.  with  levels  in  excess  of  52  UB 
filled  in  for  emphasis.  The  500  m  depth  contour  from  a  standard 
chart  is  shown  for  orientation.  «  =  0.001  dBdcm. 
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Fig.  12.  Close-in  baekscattering  levels,  experiment  14.  Contours: 
-27  to  -56  dB  in  5  dB  steps.  «  0.001  dB/km. 


or  f4  (Bass  and  Fuks).  A  review  by  Back  (1981 J  of  sound 
transmission  data  in  the  Arctic  gives  comparable  results  at  9 
Hz  (ao  =  3  x  10  ’  dB/km  for  £,  -  2  m),  but  he  finds  the  data 
to  have  a  frequency  dependence  of  /. 

Another  possibly  dominant  loss  mechanism  at  9  Hz  is  the 
third  one,  volume  scattering.  A  recent  review  [Kibbleuhiw 
and  Hampton.  1980)  has  shown  a  mean  value  aj  -  10  '  dB/ 
km  with  standard  deviation  of  the  same  order  for  seven 
experiments  in  various  regions  of  the  deep  Atlantic  and 
Pacific.  (Data  in  the  eastern  North  Pacific  for  acoustic  paths 
through  the  subarctic  convergence  were  excluded  by  the 
review  authors  as  being  anomalously  high.)  Models  exist 
which  connect  volume  scattering  loss  with  sound  speed 


fluctuations  ( Mctlcn  ct  al..  1974]  which,  in  turn,  can  he 
caused  by  internal  waves  [Mcllcn  ct  al..  1976;  Halle.  1979], 
among  amny  possible  mechanisms.  The  mean  estimate  of 
would  have  a  5  dB  effect  at  2500  km.  with  negligible  effect 
and  a  10  dB  effect  at  the  I  standard  deviation  extremes.  ct\  is 
frequency  independent. 

In  the  absence  during  CANBARX  of  direct  ice  scattering 
data,  and  data  on  internal  waves  and  other  water  column 
fluctuations  which  could  help  fix  a.  we  adopt  as  a  nominal 
value,  a  =  lO-*  dB/km.  This  value  has  the  virtue  of  being 
near  the  center  of  uncertainty,  from  which  we  might  infer 
‘better’  values  upon  inspection  of  the  baekscattering  levels. 


Fig.  II.  Close-in  baekscattering  levels  for  two  different  shots, 
both  200  kg.  Dotted  line  experiment  15.  solid  line  expeiimcnl  14. 
Contours:  28.  -  52.  V,  dB.  a  0  001  dB/km. 


5.  Backscattering  Data 

With  the  nominal  value  a  =  10  ’  dB/km.  and  other 
parameters  estimated  previously,  we  can  chart  the  backscat- 
ter  level 


BL  =  10  log  m„.  dB  (22) 

where  /a,,,  with  use  of  (21).  can  be  found  from  the  measured 
reverberation  level.  Each  of  the  backscatter  charts  we  will 
show  has  been  cleansed  of  array  side  lobes,  since  maximum 
likelihood  processing  can  introduce  spurious  values  at  angles 
other  than  those  corresponding  to  real  returns.  In  essence 
we  have  sharply  attenuated  these  side  lobes  via  an  algorithm 
related  to  the  directivity  of  strong  returns  [  Williams .  1951). 

Figure  10  shows  backscatter  from  the  Arctic  basin,  with 
strong  returns  from  some  of  the  more  distant  regions.  At  (he 
scale  shewn,  (he  chart  is  difficult  to  interpret  in  detail,  but  i( 
is  clear  (hat  CANBARX  provides  wide  area  information  on 
baekscattering. 

Contours  from  different  shots  agree  within  about  I  dB.  as 
Figure  1 1  shows.  This  figure,  for  the  Beaufort  Sea.  also 
displays  more  energetic  levels,  with  each  contour  level  4  dB 
higher  than  in  Figure  10.  Further  details  in  (he  Beaufort  Sea 
can  be  seen  in  Figure  12.  w  hich  has  the  same  lowest  contour 
as  in  Figure  II.  but  steps  up  in  finer  increments.  Between 
about  1 56 'W  and  16.5"W  the  Northwind  Escarpment  shows 


» 


> 


■V>  ’V-  ■jVwV  ‘-.v. 


8 


i‘«  >>lliV  Jl*  J  I*  j**>tl*  Ji*  A*  A^ll* J>*»*>*  .li,.i»*.i>*->t».l««.hl‘.<fcl  . 


b^Ja^Jh*  it*  i 


l)vi  r  et  at.:  Acoustic  Backscaiter  From  Arctic  Basin 


boldly,  with  Chukchi  features  showing  beyond.  To  the  south 
the  Alaska  slope  also  shows  boldly  and  continues  toward  the 
Canadian  Archipelago,  with  the  slopes  off  Tuktoyaktuk. 
Banks  Island,  and  Prince  Patrick  Island  also  showing  well 
(see  71°N,  135°W;  73CN,  I28°W;  and  77:N.  125°W.  respec¬ 
tively).  These  Beaufort  Sea  data  are  little  affected  by  the 
aforementioned  uncertainties  in  a.  since  uR  is  reasonably 
small  for  the  range  of  R  encountered. 

A  prominent  backscaltering  feature  often  omitted  in  exist¬ 
ing  topographic  charts,  is  centered  at  about  73.2°N,  139. (fW. 
While  some  charts  indicate  one  or  more  features  in  this  area 
[e.g.,  De  Lector,  1967;  National  Geographic,  1971 ;  Stefans- 
son,  1921;  U  S.  Army  el  til..  1971)  none  show  s  them  of  scale 
or  relief  sufficient  to  explain  the  observed  backscatter.  For 
example,  we  would  expect  the  feature  to  be  shallower  than 
1000  m  and/or  to  have  steep  slopes  as  in  the  Northwind 
Escarpment.  Standard  bathymetric  measurements  in  the 
area  would  provide  the  ultimate  test  of  its  reality,  but  in  the 
interim  we  are  encouraged  to  accept  it  not  only  on  the  basis 
of  the  present  result,  but  also  because  magnetic  anomaly 
data  are  supportive  [ Taylor  et  al.,  1981).  We  tentatively 
name  this  the  'G.  Leonard  Johnson  Feature'  to  honor  the 
individual  who  most  encouraged  our  work. 

In  some  ways,  data  beyond  about  1500  km  are  disappoint¬ 
ing  or  unclear  in  physical  message.  Figure  10  shows  some 
returns  on  ’land";  these  must  be  due  to  secondary  basin 
baekscattcring  paths  which  we  have  not  as  yet  unraveled. 
Inability  to  follow  the  shape  of  the  continental  slope  at  long 
ranges,  also  to  be  seen  in  Figure  10.  is  a  consequence  of 
finite  azimuthal  resolution  w  hich  is  us  gross  as  250  km  or  so 
at  large  R. 

One  may  also  question  the  transmission/back  scattering 
models  and/or  estimates  used  by  us,  at  least  because  the 
most  distant  levels  seem  only  occasionally  to  be  as  high  as 
those  close  in.  We  offer  some  discussion  points  on  this; 

1 .  A  larger  value  of  a  would  have  'brightened'  the  distant 
returns.  Enhancement  by  5  dB  or  more  at  2500  km  would 
seem  appropriate,  but  in  view  of  point  4  we  are  reluctant  to 
ascribe  all  of  it  to  a  too  low  value  of  a. 

2.  Our  model  holds  that  the  incident  energy  below  about 
1000  m  can  be  neglected.  We  see,  however,  that  the  Nort,i- 
wind  Escarpment,  believed  to  be  no  shallower  than  about 
1000  m.  comes  in  strongly.  We  interpret  this  as  evidence  of 
abnormally  large  slopes  in  this  feature,  and  therefore  our 
model  should  be  extended  to  allow  the  combination  of 
weakly  energized  but  strong  scatterers  at  depth. 

3.  The  attitude  maintained  in  point  2  on  concentration  of 
energy  above  about  1000  m,  also  helps  explain  the  returns 
from  features  beyond  the  Northwind  Escarpment  (see  vicini¬ 
ty  of  78°N,  I60°W  in  Figure  10).  That  we  have  irradiated 
them  sufficiently  to  obtain  returns  is  clear;  but  not  included 
in  the  present  transmission  model  is  loss  caused  by  the 
shielding  feature  in  the  return  direction.  A  rough  estimate  of 
the  return  loss  is  3-6  dB  over  the  Chukchi  Plateau,  and  such 
losses  applied  to  the  normalized  data  would  enhance  the 
entire  Siberian  slope,  as  well  as  particular  features  along  it. 

4.  Similar  arguments  help  explain  in  part  the  relatively 
weak  but  measurable  returns  from  the  European  Slope  and 
from  Greenland  and  the  Canadian  Archipelago  beyond  the 
Alpha  Ridge.  The  Alpha,  Lomonosov,  and  mid-Arctic  ridges 
probably  cause  return  losses  of  a  few  dB  each,  and  also  may 
time-spread  the  signals  in  such  a  way  as  to  blur  the  return 
patterns. 


6.  Speculations  on  Basin  Roughness 

In  an  attempt  to  understand  the  backscaltering  mechanism 
underlying  the  observed  tesults,  we  ascribe  to  the  slope  the 
follow  ing  characteristics:  ( I )  It  is  a  rough  surface,  statistical¬ 
ly  homogeneous  throughout  the  Arctic.  (2)  the  roughness 
height  measured  from  the  mean  surface  has  a  Gaussian  pdf 
with  £:  its  mean  square  value,  and  (3)  the  heights  have  an 
auto-correlation  function  exp  [  -v:/A:)  with  y  the  in-plane 
distance  separating  two  observation  positions  (isotropy  on 
each  wedge  surface).  We  adopt  standard  scattering  theory 
[ Beckmann  and  Spizzii  hino.  1963)  which,  in  small  angle 
approximation,  yields  the  energy  scattering  coefficient  in  the 
backscaltering  plane 


bs  =  -  A4  <<t>,  +  4),)4 

O 


■  exp - —  -  <JV>:  (23) 

16 

This  equation  is  valid  for  k:.\:  I  and  A:£':|<t>,  t  <l*/,):  «  I. 

where  k  is  (he  wave  number,  and  <J»  -  d>  +  i/>  is  the  ray  angle 
measured  from  the  wedge  for  each  of  the  incident  and 
backscatter  fields.  The  first  wave  number  condition  restricts 
attention  to  features  whose  jaggedness  is  not  abrupt  on  a 
scale  less  than  a  quarter  wavelength,  and  the  second  classi¬ 
fies  the  surface  as  only  slightly  rough,  as  all  surfaces  must  be 
near  glancing  incidence.  We  believe  these  conditions  are 
physically  sensible,  at  least  in  the  main. 

Equation  (23)  can  be  simplified  further.  The  argument  of 
the  exponential  will  be  small  if  A  A  <  4/<J>:.  From  equation 
(10)  the  maximum  value  of  d>:  is  4 ihh.  Thus  if  A  A  <  (iA0,)  1 
we  have 

bs  =  \  k4  i:\:  (ct>,  +  <L„)4  (24) 

O 

We  see  that  1  <  k.\  <  10"  would  satisfy  both  the  latter  and 
earlier  conditions  on  A  A.  The  angle  average  of  6,  is  then  easy 
to  calculate  and  is 


.  ,  ,  J  2  I  /  8,  0,. 

‘  ■  2‘ *■  T  *  5  (i  '  7, 


where  the  interval  averaged  over  is  (4 i//0,.)'':  and  (4i//0|)' :  for 
and  <J>fr,  respectively. 

Numerically  we  have,  with  0,  =  0.14  and  0,  =  0.30,  b  ~ 
(0.  l8)A4{2A2(Jr  which  we  used  in  the  Ihs  of  equation  (21). 
Thus  =  (O.WAVaV).  With  the  use  of  the  half- 
Gaussian  pdf  for  V  we  get  (i Ir')  =  np^',  and  with  the  estimate 

p»v  =  3  x  10  3  and  the  datum  pn  ~  6  x  10  4  ti  e..  BL - 32 

dB  as  a  slope  average),  we  estimate 

A;£A  -  6.3,  (26) 

or  £.\  —  4500  nr.  It  is  relevant  to  add  that  use  of  rj  directly 
from  equation  (20)  without  basin  averaging  leads  to  the  same 
estimate  of  k2C,  \  (e  g.,  p  ~  10  \  ifr  ~  4  x  I0~\  and  <Jr'  -  2  x 
10  :,  the  latter  two  at  factors  of  about  I  standard  deviation 
away  from  the  mean). 

We  can't  extract  more  than  the  estimate  of  4500  nr’  for  the 
product  of  rms  height  and  correlation  radius.  We  note. 
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however,  lhal  Ihe  rms  roughness  slope,  for  small  slopes,  is 
2l/:£/A,  and  in  ihe  range  of  ihe  model's  validity  ( 1  <  £  <  I0: 
m  and  3  x  |()'  >  A  >  30  m)  the  rms  slope  is  between  3  x  10  4 
and  some  number  Oil),  comfortably  bracketing  the  rough¬ 
ness  slope  one  might  expect. 

There  is  another  possible  insight  on  basin  backscattering. 
It  might  have  been  expected  th..<  the  rough  underside  of  the 
ice  would  be  at  least  as  impoitant  in  backscattering  as  the 
continental  margin.  But  our  data  are  dominated  by  the 
margin,  at  least  out  to  R  -  20(H)  km.  If  we  take  the  observed 
levels  in  the  center  of  the  east  Arctic  Ocean  as  an  upper  limit 
of  underice  backscattering  (  40  dB,  see  Figure  10)  we  find 
with  suitable  adjustments  that  the  ice  product  £,A,  is  no  more 
than  350  nr.  (8  dB  less  than  mean  slope  levels.  10  dB  less 
because  '.he  acoustic  window  always  fully  intersects  the  ice. 
1 1  dB  less  because  ray  angles  striking  the  ice  are  on  average 
about  1 .9  larger,  4  dB  more  because  a  is  probably  twice  the 
nominal  value.  3  dB  more  because  loss  in  round  trip  propa¬ 
gation  to  the  East  Arctic  did  not  account  for  the  Alpha, 
Lomonosov,  and  mid-Arctic  Ridges.)  We  previously  esti¬ 
mated  these  parameters  in  our  discussio  of  transmission, 
and  their  product  (170  nr’)  is  about  1/2  ot  that  deduced  here 
as  an  upper  limit.  Thus,  some  support  is  evident  for  the 
importance  of  ice  scattering,  and  for  a  value  of  a  larger  than 
the  nominal  one  used,  perhaps  as  large  as  a  -  2  x  10  '  dB/ 
km.  Also  a  possible  explanation  for  the  diffuseness  of  the 
distant  returns,  additional  to  those  offered  in  section  5,  may 
be  underice  backscattering.  Discrimination  against  such 
signals  depends  upon  and  thus  higher  resolution  systems 
could  improve  distant  returns,  to  the  extent  they  are  limited 
by  ice  scattering. 

We  do  not  wish  to  overstate  the  firmness  of  equation  (26) 
nor  the  analogous  conclusions  on  ice  backscattering.  Based 
on  absolute  error  equation  (26)  is  accurate  to  no  better  than  a 
factor  of  2(3  dB),  and  must  be  considered  significantly  worse 
given  the  guesses  we  have  had  to  make  on  i/»and  its  pdf.  The 
‘adjustments’  used  to  reach  conclusions  on  ice  backscatter¬ 
ing.  though  perhaps  reasonable,  are  pure  speculations  and 
need  more  comprehensive  research  to  gain  credence. 

Finally,  a  note  on  experiment  design  is  appropriate.  The 
expected  backscatter  is  proportional  to  f*.  so  it  might  be 
argued  that  higher  frequency  measurements  are  better,  espe¬ 
cially  given  that  array  performance  would  improve.  But 
losses  also  would  increase  and  unfortunately  might  be  as 
drastic  as  exp[-n/4],  where  n  is  some  range  dependent 
number.  Also  the  source  spectral  density  goes  as  f~2,  and 
the  noise  as  /°  to  /  2,  so  it  might  turn  out  that  higher 
frequencies  are  not  better  after  all. 

7.  Conclusions 

Compared  with  omnidirectional  systems,  an  array  can 
significantly  enrich  the  quality  of  reverberation  or  backscat¬ 
ter  data.  CANBARX  data  provide  new  insights  on  basin 
features  related  to  backscattering,  including  possible  delin¬ 
eation  of  new  or  little  explored  topographic  features,  and 
estimation  of  roughness  parameters  for  the  continental 
slope.  We  are  on  the  whole  pleased  w  ith  the  promise  of  this 
technique  and,  in  its  first  use  in  CANBARX.  have  gained 
experience  essential  to  father  development. 

Many  improvements  and  extensions  of  the  results  present¬ 
ed  here  can  be  suggested:  ( I )  Clearly,  further  analysis  of  the 
tentative  G.  Leonard  Johnson  Feature,  and  other  ‘abnormal' 
returns,  should  be  carried  out.  (2)  Backscattering  charts  for 


center  frequencies  other  than  9  Hz  should  be  prepared  to 
check  the  frequency  dependence  of  the  many  physical 
processes  involved.  (3)  More  careful  modeling  of  the  ice 
backscatter  process  would  be  useful  not  only  in  clarifying  its 
relationship  to  margin  backscattering.  but  also  for  its  intrin¬ 
sic  value.  There  is  promise  that  ice  backscattering  in  the 
manner  treated  here  can  be  used  to  gain  synoptic  informa¬ 
tion  on  ice  roughness  over  large  areas.  (4)  Statistical  analysis 
of  the  temporal  extent  of  the  signal  envelopes  can  yield  the 
pdf  of  i /»,  which  can  be  used  to  check  or  modify  the  one 
assumed,  and  thus  should  be  done.  (5)  Other  measurement 
locations  should  be  tried  to  test  the  generality  of  the  method, 
and  one  such  has  been  undertaken  [Baggeroer  and  Dyer. 
I980J.  (6)  A  scries  of  shots  at  various  depths  would  illumi¬ 
nate  other  features  in  the  basin,  and  should  be  carried  out. 
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